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Abstract. We present a detailed theoretical study of oxygen molecules in high 
rotational states (molecular superrotors) interacting with an external magnetic 
field. The system shows rich dynamics, ranging from a spin-selective splitting 
of the angular distribution over molecular alignment to an inversion of the 
rotational direction. We find that the observed magneto-rotational effects are 
due to a spin-mediated precession of the orbital angular momentum around the 
magnetic field. Analytical expressions for the precession frequency in the limits 
of weak and strong magnetic fields are derived and used to support the proposed 
mechanism. In addition, we provide the procedure for a numerical treatment of 
oxygen superrotors in an external magnetic field. 
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1. Introduction 

The control of molecular rotation by electromagnetic 
fields has proven to be a powerful tool in a number 
of rapidly growing fields of molecular science (for a 
recent review, see [1]). Aligning the molecular axis or 
the molecular angular momentum is a key requirement 
in such diverse fields as attosecond high-harmonic 
spectroscopy [2,3], photoelectron spectroscopy [4], 
controlling molecular interactions with atoms [5], 
molecules [6], and surfaces [7,8], or altering molecular 
trajectories [9-11]. 

The alignment of the molecular axis by electro¬ 
magnetic fields is a well-researched problem [12-15]. 
Polar and paramagnetic molecules can be aligned by 
the use of static electric and magnetic fields [16,17]. 
Here, the applied field interacts with the permanent 
dipole moment of the molecule and creates pendular 
states in which the molecular axis librates [18]. The 
control can be extended to molecules without a perma¬ 
nent dipole moment by the use of laser pulses [19-22] 
with intensities of the order of several TW/cm^. The 
electric field of the laser pulse polarizes the molecule 
and subsequently interacts with the induced dipole. 
Orientation of the angular momentum, i.e. the creation 
of an ensemble of unidirectionally rotating molecules, 
can be achieved by schemes like the molecular pro¬ 
peller [23,24] which uses two well-timed cross-polarized 
laser pulses, chiral pulse trains [25] consisting of period¬ 
ically applied laser pulses, or the so-called optical cen¬ 
trifuge [26,27]. In the latter technique, the molecules 
are subject to two counter-rotating circularly polarized 
fields that are linearly chirped with respect to each 
other. The resulting interaction potential creates an 
accelerated rotating trap, bringing the molecules to a 
fast-spinning state. 

Very fast spinning molecules, with rotational 
energies comparable to the molecular bond strength, 
are known as molecular superrotors. Such states can be 
created by the above-mentioned optical centrifuge. An 
alternative method is to use a train of short laser pulses 
that are separated by the rotational period [28], using 
the rotational quantum resonance effect [29]. However, 
the latter scheme is limited because of the centrifugal 
distortion of the fast spinning molecules [30]. Due 
to their high angular momenta, reaching hundreds 
of h [31], superrotors can be expected to show a 
qualitative different behaviour than normal molecular 
rotors. In pioneering work [27], it was shown that 


the centrifugal forces of the superrotors can be strong 
enough to break molecular bonds. It was also 
found [32,33] that superrotors are very robust against 
collisions, showing gyroscopic dynamics over a long 
time [34], followed by an explosion-like equilibration 
of the internal energy. 

In a recent experiment, the dynamics of param¬ 
agnetic superrotors in a static magnetic field was in¬ 
vestigated [35]. It was found that a gas of param¬ 
agnetic superrotors can become optically birefringent 
when a magnetic field is applied. This effect was called 
magneto-rotational birefringence. As an example for 
a paramagnetic superrotor molecular oxygen O 2 was 
used, which has a non-zero electronic spin of S' = 1 in 
its electronic ground state . By means of the op¬ 

tical centrifuge technique, the molecules were excited 
to unidirectional rotation with angular momenta of up 
to 100?i - ten times higher than for a typical oxygen 
molecule at room temperature. A magnetic field of 3 T, 
applied at 90 degrees to the laser propagation axis (i.e., 
the axis of rotation), lead to optical birefringence of the 
molecular gas, indicating alignment of the molecules. 
Noteworthy, the considered magnetic field strength was 
by far too low to induce the observed alignment in a 
thermal sample of oxygen molecules. As a mechanism, 
it was proposed that the magnetic field leads to a pre¬ 
cession of the large orbital angular momentum N of the 
centrifuged molecules around the magnetic field. This 
precession then leads to the alignment of the molecu¬ 
lar axis along the magnetic field. It was concluded that 
the interaction of N and the magnetic field is mediated 
via the electronic spin. 

In this article, a detailed theoretical analysis of the 
magneto-rotational effect reported in [35] is provided. 
Further insights using an experimental approach are 
provided in a companion article [36]. 

The present article is structured as follows. 
In section 2, the theoretical model is introduced. 
The effective Hamiltonian of the rovibronic ground 
state is described, as well as the numerical methods 
to calculate the rotational dynamics in the optical 
centrifuge and the probe signal. Also, the method to 
extract the angular distribution of the molecular axis 
from the spin-rotational wave function is explained. In 
section 3, the results of the numerical simulations are 
presented, in particular the simulated probe signals 
as well as the angular distributions. In section 4, 
the proposed mechanism of the magneto-rotational 
effect - spin mediated precession of the orbital angular 
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Figure 1. A sketch of the considered scenario. An optical 
centrifuge laser pulse interacts with O 2 molecules at room 
temperature. Additionally, a static magnetic field B directed 
at 90 degrees to the laser propagation axis is applied. 

momentum around the magnetic field - is explained 
in detail. An analytical expression for the precession 
frequency is derived for the limiting cases of weak and 
strong magnetic fields, and used to test the validity of 
the proposed mechanism. 

2. Model and numerical methods 

2.1. Considered scenario 

In this article we consider a scenario resembling the 
experiment described in [35]. It is depicted in figure 1: 
An optical centrifuge laser pulse [26,27] interacts with 
molecular oxygen at room temperature. As a result, 
the molecules are excited to fast unidirectional rotation 
in the polarization plane of the laser. A static magnetic 
field B is applied perpendicularly to the propagation 
direction of the laser pulse. The fast rotation of the 
molecules is measured by a second, circularly polarized 
laser pulse, which propagates in the same direction as 
the centrifuge pulse. 

In the following, we define the laboratory z-axis 
as the direction of the magnetic field, and the x-axis 
as the laser propagation direction. 

In our model, we assume that the intensities of the 
laser field and the static magnetic field are spatially 
uniform. Also, we neglect molecular collisions. 
Experimental results have shown that collisional effects 
and the magneto-rotational birefringence can be well 
distinguished [33,35], justifying this simplification of 
the model. 

2.2. Molecular Hamiltonian 

To describe the rotational dynamics of the centrifuged 
molecules, we use an effective Hamiltonian for the 


vibronic ground state. Although the optical centrifuge 
excites very high rotational states, it was found 
experimentally [31] that for the scenario considered 
in this work, the rotational levels are still very well 
described by the Dunham expansion. We therefore do 
not need to consider the vibrational states directly, but 
rather account for the non-rigidity of the molecules by 
centrifugal distortion terms. 

Molecular oxygen has an electronic ground state 
with two unpaired electrons, and thus a non-zero total 
electronic spin S; in particular, S' = 1. This gives rise 
to spin-spin as well as spin-orbit interactions. 

The effective Hamiltonian for the oxygen molecules 
in the magnetic field is given as [37-40] 

Heff = Hrot + 'Has + 'Hso + ^ms • ( 1 ) 

Here, Hrot describes the rotation of the molecule, 
Uss and Hso the spin-spin and spin-orbit interactions, 
respectively, and Hms the interaction of the magnetic 
field with the molecule, in particular with the electronic 
spin. 

The rotational Hamiltonian is given as 

Hrot = - DN^ . (2) 

The rotational and centrifugal distortion constants 
are given as Bq = 2.8400 • 10“^^ J [41] and D = 
9.612 • 10“^® J [41], where the former is for the 
vibronic ground state. Note that we use the total 
orbital angular momentum N and not the nuclear 
orbital angular momentum O. They are approximately 
equal, with small differences arising from the electronic 
orbital angular momentum [37]. 

Since there are unpaired electrons in molecular 
oxygen, the interaction between the electronic spins 
does not average to zero. Tinkham and Strandberg 
have shown [37] that the two unpaired electrons have a 
large probability to be found at the same nucleus, and 
that the spin-spin interaction is minimized when S is 
perpendicular to the molecular axis. In particular, it 
is given as 

Ha. = \\ [352 _ 52 ] ^ ( 3 ) 

where S(^ is the projection of the total electronic spin 
on the intermolecular axis and A = 3.943- 10“^^ J [40] 
is the spin-spin coupling constant. 

The spin-orbit coupling between the electronic 
spin S and the orbital angular momentum TV is mainly 
caused by the interaction of the unpaired electrons 
with the magnetic field of the electrons that follow the 
rotating nuclei [37]. It is given as 

Hao=lN-S, (4) 

where 7 = —1.674 • 10“^® J [40]. 

The by-far strongest interaction between the 
magnetic field B and the molecule is the spin-Zeeman 
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interaction between the field and the electronic spin. 
The spin-Zeeman Hamiltonian is given as 

'Hnis = gsti-sS ■ B , (5) 

with the electronic spin g-factor gs = 2.002064 [42] 
and the Bohr magneton /is. In principle, the magnetic 
field also interacts directly with the orbital angular 
momentum TV. However, the respective g-factor is 
four orders of magnitude smaller than gs and therefore 
negligible on the time-scales considered here. 

For the numerical treatment of the problem, it is 
convenient to work in Hund’s case (b) basis [40,43]. 
In Hund’s coupling case (b), the electronic spin S is 
coupled to the orbital angular momentum TV to form 
the total angular momentum J. The basis functions 
can be denoted as 

\<P) = \p,A-,N,A;N,S,J,M). ( 6 ) 

Here, ?7 is a collective quantum number for the vibronic 
state, A is the projection of the electronic orbital 
angular momentum on the molecular axis, and M 
is the projection of the total angular momentum 
on the space-fixed z-axis. Since rj, A = 0, and 
S' = I are constant, we will use the short notation 
|? 7 , A; TV, A; TV, S, J, TH) = [JTVTkf). The rotational 
Hamiltonian is diagonal in Hund’s case (b) basis, and 
the matrix elements are 

(A^JTWl'HrotjiVJTW) = BoN{N + 1)-DN^{N + lf .(7) 
The non-zero elements of the spin-spin Hamiltonian are 
{NJM\nss\N'JM) = 

^2^(_I)^+V(27v + 1){2N' + 1) 


I ^ 

TV' 

J \ 

( ^ 

2 

N' \ 

1 N 

S 

2 J 

V 0 

0 

0 J 


where TV' = TV, TV ± 2. The round brackets denote the 
Wigner 3j-symbol, the curly brackets the Wigner 6j- 
symbol. The spin-orbit Hamiltonian is diagonal again, 
with the elements 

{NJM\Hso\NJM) = 

yTV(TVTl)WTlJ yS^Sn^T)(2S^H) 

X7(-l)»««{ * " f }. (9) 

The non-zero matrix elements of the spin-Zeeman term 
are 

(NJMlHmslNJ'M) = 

J 1 J' \ ( S J' N \ 

-M 0 M ){ J S IJ 

X V(2J' + 1)(2J + 1)S(S + 1)(2S + 1), (10) 

where J' = J, J ± 1. 


The eigenfunctions and eigenenergies of the 
effective Hamiltonian (1) are calculated via 

U-i7TeffU = W. ( 11 ) 

The columns of the unitary matrix U are the 
eigenvectors \'tpn){B) in the Hund’s case (b) basis, and 
W is a diagonal matrix consisting of the eigenvalues 
En{B), where n is a collective quantum number of the 
spin-rotation states. The strength B of the magnetic 
field enters as a parameter. Equation (11) is solved 
numerically using built-in functions of the Matlab 
program package [44]. 

2.3. Dynamics in the optical centrifuge 

The optical centrifuge is created by two counter¬ 
rotating circularly polarized laser fields that are 
linearly chirped with respect to each other [26, 
27]. The result is a linearly polarized field, where 
the polarization direction rotates with an increasing 
angular frequency. The central frequency of the laser 
pulses is assumed to be far detuned from any molecular 
resonances, and the laser pulse interacts with the 
molecules via Raman-type excitations of the rotational 
levels [19-21]. The non-resonant laser pulse induces 
a dipole in the molecule via its electric polarizability 
and then interacts with the induced dipole. Averaging 
over the fast oscillations of the laser field leads to an 
effective interaction potential [45] 

V = -^E^,S^{t) [r-p{t)f 

= -UoS^{t)[r-p{t)]\ (12) 

Here, Aa is the polarizability anisotropy, Eq is the 
peak strength of the laser electric field, £{t) is the 
(dimensionless) envelope of the laser electric field, r = 
{rx,ry,rz) is the orientation of the molecular axis, p 
is the (slowly rotating) polarization axis of the laser 
electric field, and Uq = {/AaE^ /4) is the trap depth. 

The evolution of the wave function is calculated 
by solving the time-dependent Schrodinger equation, 

ih^ = {H,s + V)\4t), (13) 

using the effective spin-rotational Hamiltonian "Heff 
[(l)-(5)] and the laser interaction potential V 
from (12). It is convenient to expand using the 
eigenfunctions \'fin) and eigenvalues En of "Heff: 

(14) 

n 

We insert (12) and (14) into (13) to obtain a set 
of coupled differential equations for the expansion 
coefficients, 


( 15 ) 
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For the temporal envelope of the electric field we use 


S^{t) = 


sin 

1 

sin^ 

0 




2 ioff 


for 0 < t < ton 
for ton ^ t <C. tp toff 

for tp toff ^ t tp 

else. 


(16) 

Here, ton and toff are the turn-on and turn-off time 
of the pulse, respectively, and tp is the duration of 
the whole centrifuge pulse (including the turn-on and 
turn-off). We consider a constant intensity of the field 
(apart from turn-on and turn-off), in order to gain 
a clear picture of the magneto-rotational effect, not 
blurred by molecules lost early from the centrifuge 
trap. A realistic experimental envelope is generally 
slowly decaying [46]. The polarization vector of the 
electric field is rotating in the yz-plane, and is given as 

p(t) = ysin -tzcos , (17) 

where y and z are unit vectors along the y and z axis, 
respectively, and fi is the acceleration of the centrifuge 
rotation. Using (12), (16), and (17), the interaction 
potential can be expressed as 


V{t) = -UoS^it) 


2 ■ 2 
Cy sm 


t \ ft 

+ryr^ sm{2l3—\ + rl cos^ [fi— 


(18) 


We use the following relations to express V (t) in terms 


of the rotation matrices D 


U)*. 


( 2 )* 

00 


= ---D, 

3 3 

- ^ ^ n( 2 )* _ 

— * * J-^nn 


00 


1 

7! 

1 

7! 


D 


( 2 )* 

20 


n( 2 )* _ 
^20 


76 

1 

7! 


B 


D 


( 2 )* 

-20 

( 2 )* 

-20 


= 


-D, 


( 2 )* 

00 


(19a) 
(19b) 
(19c) 

76^'° w7=^_io (19d) 

The matrix elements of the rotation matrices in Hund’s 
case (b) basis are given as 

{JNM\D^^*\JNM) = 

X 7(2J -t 1)(2J' -b l)(2fV -b 1)(2W -b 1) 


* n( 2 )* I * n( 2 )* 

= —D\fi -I- —D^ 


J 2 J' 
—M m M' 

N' J' 1 
J N 2 


N 2 N' 
0 0 0 


( 20 ) 


Mapping (20) to the basis j^/’n) by the use of U 
[see (11)] and using (19), we obtain the matrix elements 
{f^n\V{t)\'tfm) of the laser induced potential. Finally, 
the coupled differential equations (15) are solved 


numerically to obtain |4>(t)). Note that the coefficients 
Cn{t) are constant after the turn-off of the centrifuge 
pulse, so the wave function after the centrifuge pulse is 

> h)) = E . ( 21 ) 

n 

To include thermal effects, we do ensemble 
averaging. Therefore, we solve (15) for all thermally 
populated states jV'n), and average the result weighted 
by the respective Boltzmann factors. To make the 
calculation feasible, we include only 95% of the thermal 
population. The density matrix is therefore given as 

*max 

p(t)=ArE5*l'I'*W)('I'*WI, (22) 

i=l 

where A/” is a normalization factor, |4>i(t)) is the 
solution of (13) for the initial state |4>(0)) = [f/'i), gi is 
the respective Boltzmann factor of the initial state, and 
*max is chosen such that 95% of the thermal population 
is included. 

2.4- Angular distribution 

In order to calculate the angular distribution P{9, fi) 
of the molecular axis {6 and being the polar and 
azimuthal angle, respectively), it is beneficial to use 
the Hund’s case (a) basis. In Hund’s coupling case 
(a), the total angular momentum J is formed from its 
projection D on the molecular axis and from the orbital 
angular momentum O of the nuclei. In particular, 
for an oxygen molecule in its electronic ground state, 
D is equal to the projection E of the electronic spin 
on the molecular axis, where E = 0, ±1. The 
eigenfunctions of Hund’s case (a) are the symmetric top 
wave functions [47]. The internal axis of this symmetric 
top is the molecular axis, and the projection of J on the 
molecular axis is E. The angular distribution P{0, 
is therefore obtained from the wave function in Hund’s 
case (a) by integrating over the internal coordinate, 
which yields 

P{0,4’)= E 

X e},M's(70jMs(7<i>7'(7<f>M(</>), ( 23 ) 

where m'y: elements of the density 

matrix in Hund’s case (a) basis and 



(24a) 

$m(</>)= 

(24b) 


where (^) reduced Wigner matrix elements. 

They are given as [40] 
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^ V(J + M)!(J-M)!(J + S])!(J-S])! 

^ (J + M-s)!(J-S-s)!s!(s + E-M)! 

cos-j (^^^ 2 ) ■ 


The index s runs over all values for which the 
arguments of the factorials are not negative. 

Unfortunately, (25) is numerically hard to evalu¬ 
ate for most values of M due to the large factorials 
involved. Therefore, the following recursion relation is 
used: 




E- (M- l)cos6» 
sm0y/{J - M +1){J + M) 




^(j + M-1)(J-M + 2) (j) 




- 2 - 


{J-M + 1){J + M) M-2.SV 

(26a) 


E- (M-h l)cos6» 


sin6lv^(J-hM-h 1)(J- M) 


' (j_M-1)(J + M + 2) 

(J-f M-hl)(J-M) )■ 

(26b) 


This recursion relation is only stable for [48] 

^ , 2MEcos6»-M2-E2 ^ 

J{J + 1) +--> 0 . 

Sin 0 


(27) 


When the inequality (27) is not fulfilled, the recursion 
has to be done in a direction which leads towards 
fulfilling the inequality [48]. For E = 0, this means 
using the recursion relation (26) with decreasing |M|, 
in particular (26a) for negative M and (26b) for 
positive M. This scheme can also be used for the other 
two cases of E = ±1 with an error negligible for our 
purposes. The starting points for the recursion are 
[obtained from (25)] 


4 +nE('^) = '^-]-i.i:(^) = 0 (28a) 



(28c) 


A test of the accuracy of the numerical method via 
the orthonormality of the obtained symmetric top 
eigenfunctions showed an error of less than 100 ppm, 
for angular momenta up to J = 42. This is sufficient 
for our purposes. 

To obtain the angular distribution from (23), one 
needs to map the density matrix (22) to the basis of 


Hund’s case (a). In particular, the mapping is given 
as [40] 

(jMEijfvM) = {-ly-^yw+i ^ ^ ^. 

(29) 


2.5. Raman signal 


The ultrafast rotation of the centrifuged molecules 
can be directly measured by Raman spectroscopy, 
using a circularly polarized probe pulse [31]. The 
coherence between states with AN = ±2 causes a 
Raman shift of twice the rotational frequency of the 
centrifuged molecules. The shift is negative (positive) 
if the molecules rotate in the same direction (opposite 
direction) as the probe pulse. This frequency shift can 
be interpreted as a rotational Doppler effect [49]. 

The probe signals are calculated using the model 
described in [50]. The output Flout of the probe field 
(with left and right circular polarized components) is 
given as 


Eoutid) « e*“^*exp(zAGd) 


X 


cos {I3i\p\d)l -t-zsin(/3i|p|d) 


0 e-*'*' 

0 


Ein, 

(30) 


where Flin is the input of the probe field, d is the 
length of the gas cell, fit is a parameter depending 
on the experimental conditions, LOi is the central 
frequency of the probe, F is the identity matrix, G 
is the expectation value of the alignment to the laser 
propagation axis, and p = \p\e^^ can be understood 
as a “circular alignment” in the centrifuge plane. We 
use dfii = 10 in our calculations, as this approximately 
matches typical experimental conditions [31]. For the 
coordinate system from figure 1, the quantities G and 


p are given as [50] 

G = {rl) (31a) 

p ={{ry+irfi)y. (31b) 

The simulated probe signal is obtained from (30). The 
necessary expectation values are obtained as 

G{t) =Tr{Gg{t)} (32a) 

pit) = Tr{pp(t)} . (32b) 


Here, g is the density matrix, and G and p are the 
matrix representations of G and p [see (31)]. The 
matrix elements of G and p are obtained by the help 
of (31), (19) and (20). As input probe we use a 
clockwise polarized pulse with a sine square envelope. 
The measured signal is the anticlockwise part of the 
output of the probe. It is given as a function of the 
pump-probe delay tp as 

S{t, tp) = sin (10|/9|) - tp) , (33) 
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where fin is the envelope of the probe pulse. Finally, we 
take the Fourier transform over t to obtain the output 
spectrum |f(w,tp)p. 

3. Results of the simulations 

In this section, the results of the simulations are 
presented. The parameters were chosen such that 
they closely resemble the experiment described in [35]. 
We consider a centrifuge pulse that is anticlockwise 
rotating around the x-axis, and therefore excites 
a rotational wavepacket with the orbital angular 
momentum N pointing along —x. The duration of 
the laser pulses is tp = 36 ps, with an acceleration of 
/3 = 0.52 • 10^'^ s“^ and a turn-on and turn-off time of 
ton = toff = 2 ps. The intensity of the laser pulse is 
Iq = 2 TW/cm^, which corresponds to a trap depth 
of Uo ~ 16 Bo (where Bq is the rotational constant). 
At the end of the pulse, the centrifuge is spinning 
with 3 THz, which corresponds to a rotational level 
of Ri 35. The initial rotational temperature of the 
molecules is set to room temperature, T = 298 K. 
The strength of the magnetic field is R = 3 T, and 
its direction - defining the z-axis - is perpendicular to 
the laser-propagation axis. The simulated probe pulses 
are right-circularly polarized, and have a duration of 
10 ps (full width at half maximum), corresponding to 
a resolution of approximately 3 cm“^. 

Figure 2 shows the simulated spectrum. For 
comparison, also the case of no magnetic held 
is presented. The simulation reproduces recent 
experimental results very well [31,33,35]. During 
the centrifuge pulse (the hrst 36 ps), whilst the 
molecules are accelerated, one can see a strong anti- 
Stokes signal (the centrifuge and the probe pulse 
are rotating in opposite directions) with a linearly 
increasing frequency shift. After the end of the laser 
pulse, the signal remains at a constant shift and 
consists of two lines around 200 cm“^. This shift 
is caused by molecules rotating anticlockwise around 
the x-axis. In particular, the two lines correspond to 
the levels N = 35 and N — 33 (remember that the 
shift is twice the rotational frequency, see section 2.5). 
After about 500 ps, the Raman signal fades. This 
is due to the coherent dephasing of the three spin 
states [33]. For the case of no magnetic field, the 
Raman signal partially revives after about 1 ns, and 
a full revival can be observed after about 2 ns (not 
shown in figure 2). If a magnetic field is applied, no 
such revivals are observed. However, another signal 
appears: After about 300 ps, a Rayleigh component 
becomes visible. Such a line has also been observed in 
a recent experiment [35]. 

Figure 3 shows the strength of the three parts 
(Stokes, Rayleigh, and anti-Stokes) of the probe signal. 



100 200 300 400 500 600 



Figure 2. Probe spectrum for centrifuged oxygen molecules. In 
panel (a) the case of no magnetic field is shown, in panel (b) the 
case for a magnetic field of 3 T. The colour axis is in arbitrary 
units. See text for the parameters. 



pump-probe delay in ps 


Figure 3. Relative strengths of the probe signal components 
(anti-Stokes, Rayleigh, and Stokes) for centrifuged molecules 
in a magnetic field of 3T. The strong peak of the anti-Stokes 
component during the centrifuge pulse (the first 36 ps) is due 
to the strong alignment in the centrifuge field. See text for the 
parameters. 
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for pump-probe delays of up to 2 ns. It can be seen that 
the anti-Stokes signal decreases almost monotonically, 
with a weak local maximum at around 1300 ps. 
The Rayleigh signal slowly builds up, reaching a 
maximum after 400 ps, then remaining approximately 
constant. Remarkably, also a Stokes signal builds up 
over time. It reaches a maximum at 800 ps and is about 
hundred times weaker than the other two signals. The 
appearance of the Stokes signal indicates that some 
of the molecules may have changed their direction of 
rotation. An experimental observation of this Stokes 
signal has also been reported in [35]. For the case of no 
magnetic field (not shown here), neither a Stokes nor 
a Rayleigh signal can be observed. 

Additional insights into the dynamics can be 
gained from the angular distribution. It is shown 
in figure 4 for different times after the end of the 
centrifuge laser pulse (a short movie can be found in 
the supplementary material). For no magnetic field, we 
can see that the distribution is elongated along an axis 
that rotates anticlockwise around the laser propagation 
axis. The dynamics does not change significantly over 
time [compare figures 4 (a) and (b)]. The picture 
is quite different when a magnetic field is applied. 
Shortly after the centrifuge pulse [figure 4 (c)], we 
can see approximately the same distribution as for 
B = 0. At later times [figure 4 (d)], a very different 
picture emerges. Whilst the angular distribution is still 
elongated, it now points along the magnetic field and is 
not rotating any more. This alignment along the z-axis 
can explain the appearance of the Rayleigh component 
in the probe signal. 

4. Rotation of the angular momentum by the 
magnetic field 

In order to explain the observed birefringence and 
Stokes signals, it was proposed in [35] that the orbital 
angular momentum TV precesses around the magnetic 
field B. This scenario is depicted in figure 5 (a). 
Right after the centrifuge pulse (leftmost panel), N 
points along —x and the molecules rotate in the yz- 
plane. A magnetic field induced torque r acts on TV, 
leading to a precession of TV around B. The situation 
as seen for the probe pulse is shown in figure 5 (b), 
where the projection onto the yz-plane is shown. At 
the beginning (leftmost panel) a anticlockwise rotation 
of the molecules in the j/z-plane is seen. Accordingly, 
the probe signal has a strong anti-Stokes component 
(remember that the probe is right-circularly polarized). 
After some time, when TV points along the y-axis, the 
molecules rotate in the a;z-plane. Projected onto the 
yz-plane, this is seen as alignment along the z-axis 
(middle panel). Accordingly, the probe signal acquires 
a Rayleigh component. Eventually, TV points along 


+x and uni-directional rotation is seen again in the yz- 
plane - but now in a clockwise sense (rightmost panel). 
This leads to the observation of a Stokes component in 
the probe spectrum. 

The torque t can not originate from a direct 
Zeeeman interaction between B and TV. The 
respective g-factor is of the order of 10“^ [42]; for the 
considered field strengths of a few Tesla, this would 
correspond to precession periods of 100 ns and longer, 
clearly too slow to explain the observed effect (such 
a slow precession has been observed experimentally 
for nitrogen superrotors [36]). Instead, the interaction 
between TV and B is mediated by the electronic spin 
S. The spin interacts with the magnetic field via the 
Zeeman interaction 'Hms 7 and with TV via the spin- 
rotation interactions "Hss and "Hso- As a result, the 
effective Hamiltonian depends on the orientation of TV 
with respect to B, giving rise to the torque r. 

In order to test the proposed model, we calculate 
the frequency of the precession of TV around B. The 
torque r acting on TV is given as 


T = —TV X (Vat • WN,a{B)) , 


(34) 


where is the potential energy arising from the 

magnetic field, for a molecule in the rotational level TV 
and spin-state cr (We introduced the new label a for 
the spin-state, as the meaning of a changes with the 
strength of the magnetic field). Since the magnetic 
field points along the z-axis, we obtain from (34) the 
precession frequency h'N,a as 


= T' 
h 


1 OWnAB) 


ON, 


(35) 


where h is Planck’s constant, and TV^ is the projection 
of TV on the z-axis. In the following, we will derive 
analytical expressions of as a function of the field 
strength and the terminal velocity of the centrifuge, for 
two limiting cases: The case of a weak magnetic field, 
i.e. when the spin-Zeeman interaction is mnch weaker 
than the spin-rotation interaction, and the opposite 
case of a strong magnetic field. 


4-1. Weak magnetic field 


For a weak magnetic field, the electronic spin is 
coupled strongly to the orbital angular momentum; the 
magnetic field acts as a perturbation to the precession 
of S around TV. The spin-rotation interaction is 
therefore independent of the orientation of TV and the 
only non-zero term in i'n,<t arises from the spin-Zeeman 
term. 


gsf-B d{B- S) _ gsp.B p dS^ 
h dN^ ~ h dNA 


(36) 


Here, is the projection of S on the space-fixed z- 
axis. Averaging over the fast rotation of S around TV, 
we obtain = SnNz/\^\, where Sn ~ 0, ±1 is the 
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(a) early time, 
no B-field 







(b) early time, 
with B-field 





(c) late time, 
no B-field 




(d) late time, 
with B-field 







Figure 4. Angular distribution of the centrifuged oxygen molecules at different times, (a,b) without magnetic field and (c,d) with 
a field of S = 3 T. The red arrow depicts the propagation axis of the centrifuge pulse, the blue arrow the magnetic field. Time 
evolves from top to bottom rows, with time steps of 25 fs per row. The first row is at t = 36 ps (right after the end of the centrifuge 
laser pulse) for early times (a,c), and at t = 1000 ps for late times (b,d). See text for the parameters. See also the movie in the 
supplementary material. 































Figure 5. (a) A sketch of the proposed mechanism of the magneto-rotational effect. The panels show the situation at different 
times (time progressing from left to right). Shortly after the end of the centrifuge laser pulse, the orbital angular momentum N 
points along —x, and the molecules rotate anticlockwise in the yz-plane. A magnetic field induced torque r causes a precession of 
iV around the magnetic field B. The plane of rotation of the molecules (indicated by the blue ring) precesses accordingly. In panel 
(b), the projections of the pictures of panel (a) on the y 2 :-plane (the polarization plane of the probe pulse) are shown. It can be seen 
that at different times the molecules rotate clockwise or anticlockwise, or they appear to be aligned along the magnetic field. 


projection of 5 on TV (note that Sm characterizes the 
three spin states, thus a = S^). Thus, the precession 
frequency of TV around the magnetic field is given as 


~ R 


(37) 


Equation (37) shows that the precession frequency 
is proportional to the ratio B/N. This stems from 
the fact that the torque, induced by the spin-Zeeman 
interaction, is proportional to the strength of the 
magnetic field, but independent of TV; yet, a higher 
angular momentum needs a stronger torque to process 
with the same frequency. It can furthermore be 
seen that the torque is proportional to Sn, thus the 
three spin-states process differently: The J = TV ± 1 
states process with the same period, but in opposite 
directions, and the J = TV state does not process at all. 
Also, I’M,a is independent of TV^, thus the distribution 
of the angular momentum for a single spin state does 
not disperse, but remains in a narrow cone that was 
created by the centrifuge. 

For the angular distribution, averaged over the 
fast rotation of the molecular axis, we can conclude 
the following: Right after the centrifuge pulse, it has 
the shape of a disc in the yz-plane. In the time- 
scale of hundreds of picoseconds, this disc splits into 
three discs, corresponding to the three spin-states: One 
remaining in the yz-plane, and two rotating in opposite 
directions around z. In a very recent experiment, 
Milner et al. [36] succeeded to observe this behaviour. 

Since there is no dispersion, the Rayleigh signal 
should show clear beating with a frequency of 2vM,cr 
(only the alignment, not the orientation of TV is 
important for the signal). The signal should first peak 



magnetic field in T 


Figure 6. The time Tr of the first peak of the Rayleigh signal 
as a function of the strength B of the magnetic field, for weak 
field strengths. The solid line shows the predictions by (37) (with 
N = 35), the circles are the results of simulations using the same 
parameters (apart from B) as in section 3. 


at Tr = 1 / {4:h'N,a) after the end of the centrifuge pulse, 
when TV points along the y-axis. We compare Tr as 
predicted by the analytical formula (37) to numerical 
simulations via the Schrodinger equation (13). This 
comparison is shown in figure 6, where the solid line is 
the analytical prediction and the circles are the results 
of the simulations. A very good agreement can be 
seen. The simulations also confirmed the beating of 
the Rayleigh signal (not shown here). 

If the model is correct, we should expect a 
similar behaviour for the Stokes signal, with beating 
at VN,a and a first peak at Ts = l/( 2 i/ 7 v,cr)- 

However, simulations revealed that this is not the case. 
Furthermore, the signal is two orders of magnitude 
smaller than the anti-Stokes signal. The likely reason 
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is that the Stokes signal is - differently to the Rayleigh 
signal - affected by coherent dephasing of the different 
spin states. The Stokes signal might also be mnch more 
sensitive to weak perturbations by the magnetic field 
that are not covered by the approximations we made 
in deriving (37). 


4-2. Strong magnetic field 


For a strong magnetic field, the spin is strongly 
coupled to the field and precesses around it. The 
three spin-states are defined by the projection of S 
on the magnetic field axis, Sz = 0, ±1 = a. The 
orbital angular momentum is indirectly affected by the 
magnetic field via the spin-orbit coupling Hso and the 
spin-spin coupling Hss- The latter includes the orbital 
angular momentum implicitly via the projection S(^ of 
the spin on the molecular axis. Averaging over the fast 
rotation of the molecular axis yields S(; ~ O.Ssin^a, 
where a is the angle between S and the nuclear orbital 
angular momentum O. Approximating O hy N [37], 
we can express the spin-spin interaction as 


-Hss = xA [3A^" - 52 


1 — cos^ a — -S^ 

(j 


1 - 


(5-iV)2 _2 2 
3“^ 


(38) 


N{N + 1) 

Adding the spin-orbit coupling, we obtain the spin- 
rotation coupling between S and N as 

(S-N)^ 2 


n 


yN-S+X 


1 _ ^ _ 

^(^ + 1 ) 3 


(39) 


Averaging over the precession of S around the 
magnetic field B (which is much faster than the 
precession of N around B) yields 


n 


sr 


yNzSz - A 


TV? 


N{N + 1) 





(40) 


Since Sz and therefore the spin-Zeeman term is 
constant, (35) becomes 




I dTisr ^2 ,X Nz /3 2 _ 

hdNz h ^ hNiN + 1) V 2 ^ 


.(41) 


The precession frequency of TV for the limit of a 
strong magnetic field includes two very different terms. 
Both terms are of the same order of magnitude for 
typical oxygen superrotors (30 ^ TV < 100). The 
first one, stemming from the spin-orbit coupling, looks 
similar to the weak-field case. It is again proportional 
to the approximate quantum number defining the spin 
state. There is no 1/TV dependence, since the spin- 
orbit coupling is proportional to the size of TV, and 



500 1000 1500 2000 

pump-probe delay in ps 


Figure 7. Shown are the relative strengths of the (a) Rayleigh 
and (b) Stokes signal for centrifuged molecules, for the case of 
strong magnetic fields (in particular, B = 6,9,12 T). Note 
that the three shown Rayleigh signals are almost identical. The 
parameters (apart from B) are the same as in section 3. 


thus the torque increases with TV. The second term 
in (41) stems from the spin-spin interaction. The 
latter depends quadratically on Nz, and thus the 
interaction does not average to zero for the 52 = 0 
state. Therefore, in the strong field case all three 
spin states show a precession of TV. Additionally, the 
dependence of VN,t 7 on Nz means that the distribution 
of TV disperses during the precession, leading to a 
more or less equal distribution of TV in the xy-plane 
over time. There is no dependence on the strength 
of the magnetic field in (41). This is not surprising, 
since the torque arises from the spin-rotation coupling; 
this coupling is caused by intramolecular interactions 
and therefore independent of any external field. An 
experimental observation of the independence of the 
magneto-rotational effect from B for strong fields was 
reported in [35]: There, the birefringence signal was 
measured 1 ns after the centrifuge pulse as a function 
of the magnetic field strength. It was found that the 
signal strength becomes independent of the magnetic 
field strength for sufficiently strong fields. 

We again test the predictions of (41) by comparing 
them to simulations. The results for the Rayleigh 
signal are shown in figure 7 (a), for the same 
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parameters as in section 3, except for a stronger 
magnetic field. Changing the field strength has almost 
no influence on the signal. It furthermore remains more 
or less constant 500 ps after the centrifuge pulse, which 
is in accordance with N being equally distributed 
in the xy-plane. There is a strong peak at 400 ps 
which we can not explain with the model. The Stokes 
signal, shown in figure 7 (b), is also approximately 
independent of B. 

5. Discussion and Conclusion 

In [35] it was reported that paramagnetic superrotors 
interacting with a magnetic field B become aligned to 
the field. As a mechanism, spin-mediated precession of 
the orbital angular momentum N (which is mainly the 
angular momentum of the nuclei) around the magnetic 
held was proposed. In particular, the electronic spin 
S interacts with the magnetic held B via the spin- 
Zeeman interaction, and with the orbital angular 
momentum TV via the spin-rotation coupling. Thus, a 
strong torque stemming from the magnetic held acts on 
the orbital angular momentum. This torque is orders 
of magnitude larger than the torque stemming from a 
direct Zeeman interaction between TV and B. 

In this article, we conducted a thorough theoreti¬ 
cal analysis of the magneto-rotational effect. Starting 
from the proposed model, we derived analytical expres¬ 
sions for the precession frequency in the limits of weak 
and strong magnetic helds, respectively. This allowed 
us also to make predictions about the distribution of 
TV and the molecular axis. For weak helds, a split¬ 
ting of the angular momentum distribution into three 
parts (one for each spin state) was predicted, and the 
precession frequency is proportional to the strength of 
the magnetic held. For strong helds, it was shown that 
the angular momentum becomes evenly distributed in 
the plane perpendicular to the magnetic held, and the 
precession is independent of the held strength. A com¬ 
parison of the analytical derivations on one hand, and 
numerical simulations (reported here) as well as exper¬ 
imental observations (reported in [35,36]) on the other 
hand, showed a quantitative agreement for the weak 
held case and a qualitative agreement for the strong 
held case. 

The control of molecular rotation by electromag¬ 
netic helds is used in a large number of applications [I]. 
The magneto-rotational effect of paramagnetic super¬ 
rotors adds an additional control knob for shaping 
the angular momentum distribution of molecules. In 
particular for weak helds, a high degree of control is 
possible, including selective control of different spin 
states. The analysis of the magneto-rotational ef¬ 
fect provided in this article allows for proper tailoring 
of magnetic and centrifuge helds in possible applica¬ 


tions. The magneto-rotational effect shows qualitative 
changes for different held strengths; this makes the con¬ 
trol of molecular trajectories by inhomogeneous mag¬ 
netic helds [9-11] a particularly interesting application 
for the magneto-rotational effect. 
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